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Abstract. 

We derive the kinetic equation for pure gluon QCD plasma in a general way, 
applying the background field method. We show that the quantum kinetic equation 
contains a term as in the classical case, that describes a color charge precession of 
partons moving in the gauge field. We emphasize that this new term is necessary for 
the gauge covariance of the resulting equation. 

Heavy ion collisions at ultra-relativistic energies are widely expected to be a 
laboratory to study the formation and properties of highly excited QCD matter in 
the deconhned Quark-Gluon Plasma (QGP) phase |ij. The QGP is considered as a 
partonic system being at (or close to) local thermal equilibrium. Thus, to study the 
conditions for the possible formation of QGP in heavy ion collisions one needs to address 
the question of thermalization of the initially produced partonic medium |2J. On the 
theoretical level this requires the formulation of kinetic equations f3| involving color 
degrees of freedom and the non-Abelian structure of QCD dynamics. Different models 
for the initial conditions in ultra-relativistic heavy ion collisions suggest that at the early 
stage the medium is dominated by gluon degrees of freedom J3J; the kinetic equation for 
the pure gluon plasma is thus of special interest. 

The kinetics of a particle moving in a non-Abelian gauge field was first studied for 
SU(2) symmetry [|]. The equation of motion derived in || made it clear that in the 
presence of the gauge field the color charge of the particle keeps rotating or undergoing 
precession with time. Consequently, there is an explicit term in the classical kinetic 
equation, which describes the time evolution of the color charge of partons. Applying a 
typical method to derive a semi-classical transport equation from the quantum theory, 
which is contained in the formulation of the equation of motion for the Wigner function 
and its subsequent gradient expansion 10-0, one loses the explicit presence of the color 
charge precession because the obtained transport equations are not most simplified. The 
author of Ref. [9| simplified the kinetic equation for gluons based on Ref. [[|. However 
the simplification is under the assumption of a specific form of the covariant Green 
function (actually it is a special case of this work), and the color precession term was 
not explictly and transparently interpreted. 
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The usual treatment of the gluon transport equation is based on the decomposition 
of the gluon held into a mean held and a quantum fluctuation, or into a soft and a 
hard part. Under this approximation the gluon transport equation then describes the 
kinetics of the quanta in the classical mean held BE3- This picture is somehow similar 
to the one used while studying the energy loss of a fast parton moving in the soft mean 
held [|TT], [T5J . Since this decomposition is not rigorous, the explicit gauge invariance 
of the formalism remains unclear. To include the classical chromoheld into QCD in a 
proper way, one uses the background held method of QCD (BG-QCD) introduced by 
DeWitt and’t Hooft 11,111 in. The advantage of BG-QCD is that it is formulated in 
an explicit gauge invariant manner. 

One of the hrst attempts to derive the gluon transport equation in BG-QCD was 
presented in |l(J. But the obtained equation is not transparent. The most recent work 
has been done by Blaizot and Iancu |[1 [18| in the context of the close-time-path (CTP) 
formalism. There, however, the authors focus on formulating the transport equation in 
the vicinity of equilibrium. 

In this work, we use the CTP method to derive the kinetic equation for the gluon 
plasma in BG-QCD. Our derivation is going beyond previous results as it is quite general 
and not limited to the vicinity to equilibrium. One of the most important features of the 
kinetic equation derived here is that a color charge precession term is explicitly present, 
similar to the one in the classical equation 13 § 0 . We demonstrate that this term is 
necessary to guarantee the gauge covariance of the resulting Vlasov equation. 

In the following we use g, JiV = diag(l, —1, —1, —1) as the metric tensor, and we 
write Lorentz indices as subscripts and color ones as superscripts for the relevant 
quantities. For the gauge field and its field tensor we also denote: A M = A“T a and 
F /t „ [A] = Ff u [A]T a , where (T a ) lj = if mj are the generators of the SU( 3) c adjoint 
representation. The two-point Green function (GF) or self-energy (SE) are treated as 
matrices, and thus the color and/or Lorentz indices are sometimes omitted. 

Applying the background field method, we decompose the conventional gluon field 
into the sum of a classical background part A and a quantum fluctuation 0. Including 
the appropriate gauge fixing and ghost term for the background gauge D^[A]Qf = 0, 
the BG-QCD Lagrangian reads |15| 


C = - \f;aa + Q]F;„{A + q\ - ±(d‘;[a]q;) 2 


C‘D’M]D*{A + Q]C k 


(i) 


where D l j[A(x)] = [d x ^ — igA^{x)} 1 ^ is the covariant derivative, C l /TT are the ghost/anti 
ghost fields and a is the gauge fixing parameter. 

The above Lagrangian is invariant under the local gauge transformation of type 
I (type II transformation is irrelevant to the current problem) where the background 
field transforms as a conventional gauge fields, A/ = UA^U -1 + ig~ lTT ^ 71-1 
the gluon and ghost held transform as a matter held, Q' = UQ^U~ 

U(x) = exp [iguj a (x) T a ] is the transformation matrix. 


r_1 lira 
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The non-equilibrium dynamics is usually described in the CTP formalism |T9 


Here the action can be written as Sctp = S + — S_ + K(A±, Q±) where all fields in 
S± are defined on the positive/negative time branches and K(A±,Q±) is the kernel 
incorporating initial state correlations. The GF for the gluon has four components: 
[G ++ , G +_ , G~ + , G ] = [G F , G < , G > , G F ], marked by the positive or the negative time, 
respectively. At the tree level there are only vertices with all positive or all negative 
time variables, which differ by a sign. The non-local kernel K(A±,Q±) appearing in 
Sctp is non-vanishing only at the initial time and plays the role of the initial density 
matrix. If the initial time is not in the remote past, it introduces a time dependence 
to the distribution function, and consequently the time-translational invariance is lost 
2Cf| . It is interesting to note that the situation here is similar to the pinch singularity 


discussed in [21 


In the CTP formalism, choosing the physical representation for the GF and SE 
the Dyson-Schwinger equation (DSE) for G c reads 


D{x\)G c (aq, x 2 ) = — J d 4 x l [H c (xi,x , )G A (x , ,x 2 ) 

+ n R (aq, x')G c (x', x 2 )\ , (2) 

G c (x i,X 2 )T>*(x 2 ) = - J d 4 x'[G R (x 1 ,x')U c (x',x 2 ) 

+ G c (x 1 ,x')U A (x',x 2 )} , (3) 


where the differential operators T> and ZF in the Feynman gauge [§] (cc = 1) are expressed 
by: 

= g pa D h ;[A\D;[A] + 2 gf hai F a pcj [A\ , (4) 

®S,“ = g^D^lAlDjrlA] + 2 gf M F;,[A] , (5) 

where Djf J [A(x)\ = [d X n +igA fJj (x)Y : ’ is the conjugate covariant derivative where the 
differential operator acts on the function on its left. 

In the evolution of the gluon held one distinguishes different scales, which 
characterize quantum and soft collective motion. We introduce a mass parameter, p, as 
the separation point of the quantum and the kinetic scale. In the weak coupling limit 
g ■ Cl, the scale of collectivity ~ 1 /(g/x) is much larger than the typical extension of 
hard quantum fluctuations ~ l//i. The effect of the classical held A on the hard quanta 
involves the coupling gA to the hard propagator and is of the size of the soft wavelength 
~ 1 /(gn). The above separation of scales is the basis for the gradient expansion where 
one expresses all 2-point GF’s in terms of the relative y and the central X coordinate. 
Here are some typical scales: y — x± — x 2 ~ 1 //x, X = (aq + x 2 )/2 ~ l/(gy), A(X) ~ /i 
and F[A(X)] ~ gy 2 . 

f In the physical representation the elements of GF and SE matrices are expressed by the symmetric 
(C), the retarded (R) and the advanced (A) components. They are related to the original components 
by the following transformations: G A = G F — G > , G R = G F — G < , G c = G > + G < (the same relations 
hold for the SE IF / R / A up to a sign in front of IF and IF). 

§ In fact due to gauge conditions: [A(xi)]Gj]f(xi,x 2 ) = 0 and GJf„(xi, x 2 )Dl’ : > k [A(x 2 )] = 0, Eqs. 

( 2|,|j) are independent of the gauge parameter a. 
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In order to obtain the gauge-covariant kinetic equation one uses the gauge-covariant 
Wigner function G(X, y ) (G c , G > or G K ) defined by 

G(x u x 2 ) = V( Xl ,X)G(X, y)V(X, x 2 ) , (6) 

where V(zi,z 2 ) = T P exp(ig Jp. Z2 dz^A^) denotes a Wilson link with respect to the 
classical background held. One can also define the Wilson link as a functional of A + Q, 
but it is a much more complicated case and beyond the scope and subject of this work. 

The covariant Wigner function G(X,y) transforms as U(X)G(X,y)U~ 1 (X)] only 
the transformation at a single point X is thus relevant. For G(xi,x 2 ), however, the 
gauge transformation involves two points and therefore is not gauge-covariant. 

Taking the difference between Eq.(^j) and Eq. ([][) , we derive the following kinetic 
equation in terms of the gauge covariant Wigner function G c (X,y): 

q ■ d\G^ + ig(G^q ■ A — q ■ AG ^ 7 ) 

+ ^gquF uX (d qX G^) + ]^g(d qX G c ai )q v F uX 

+ g(F a pG - GapFfo) = 0 , ( 7 ) 

where G^ = G^yiX, q) is the Fourier transform of G(X, y) with respect to y: G(X, q) = 
f d 4 y G(X, y)e iqy . 

The above equation is located at the collective coordinate X and is gauge-covariant 
under the local gauge transformation U(X), i.e. it transforms as U{- ■ -)U~ l . Indeed, 
noting that both F /t „ and G are gauge-covariant and d q/l does not affect U, it is obvious 
that the last three terms are gauge-covariant. With some algebra, we can also prove the 
first two terms preserve gauge covariance too. 

The quantum kinetic equation Eq.(7[) is derived in a quite general manner without 
requirements or assumptions that the QCD plasma appears close to equilibrium. In 
the following we show that Eq. (0) is a natural quantum generalization of the classical 
Vlasov equation. In particular the color charge precession will be explicitly identified in 
the quantum description of color charge kinetics given by Eq.(|7j). 

The classical kinetic equation for the color singlet distribution function f(x,p,Q ) 
has the following form B 1,0: 

P - gQ a F^d pu 

- gf abc A b fl (x)Q c d Q a]f(x,p : Q) = 0 , (8) 

where Q a is the classical color charge and a = 1,.., Nf — 1. 

Comparing Eq.(0) with the quantum expression Eq.(0) it is clear that the color 
singlet distribution function / is replaced by the gauge-covariant Wigner function G c 
which is a color matrix in the adjoint representation. One can also recognize that the 
Erst, third and fourth terms of Eq.(0) are the quantum generalization of the first two 
terms in Eq.@. The last term in Eq.([7]) appears from the covariant operators [J] and 
hence is not present in the classical equation. 

|j This term can be written in different form by using generators of the Lorentz transformation in vector 
representation. Similar term can be found for the quark, but expressed through generators in spinor 
representation. 
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Particularly interesting is the appearance of the second term in Eq.(j7i). We have 
seen that its presence is crucial to assure the gauge covariance of the Vlasov equation. 
In addition this term has an interesting physical meaning. It is the quantum analogue to 
the color charge precession in the classical kinetic equation. To see this more clearly, one 
can expand G^JX,q) with respect to the generators T a of the adjoint representation 
as follows 

G c ap (X,q) = N al3 (X,q) + T‘N‘„(X,q) 

+ T‘T b N‘ b (X,q) + --- , (9) 

where T a are quantum analogues to the classical color charges Q a ; N a g(X,q ) is the 
color singlet and N“p(X,q) the color octet function, etc. This expansion is similar to 
the multipole expansion of an electromagnetic source in electrodynamics. In the lowest 
order approximation, one can keep only the Erst two terms in Eq. ([]) and neglect all 
others. Then the second term of Eq.(|7]) becomes 

!9(GV -A-q- AG %) ~ -gf^q^rdr-G^ (10) 

which is just the classical color precession term in Eq.(||). 

Since we know that D x ~ gA(X) ~ g/i , the first two terms of Eq.(|7]), i.e. q ■ d x 
and the color precession term, are at leading order 0(gg 2 ), while other terms are at 
subleading order 0(g 2 g 2 ). In the vicinity of equilibrium the natural scale in the system 
is the temperature T. The mean distance between particles is of the order of ~ 1/T, 
while 1 /(gT) characterizes the scale of collective excitations [0.0. For small coupling 
constant g these two scales are well separated. The covariant Wigner functions can 
be expanded around their equilibrium values: G = G ^ + 5G, where the equilibrium 
function G ® is a color singlet and the fluctuation 5G ~ g 2 G. Typical scales are q ~ T, 
D x ~ g 2 T , gF ~ (D x ) 2 ~ g 4 T 2 . Thus, at leading order, only the first term of Eq.flT]) 
survives and the precession term vanishes due to the color-singlet nature of G^f 

In summary, by applying the closed-time-path formalism we have derived the gluon 
kinetic equation in the background gauge of QCD. Equation (|7|) is quite general and 
is not limited to the vicinity of equilibrium. It is a kinetic equation with respect to a 
gauge-covariant Wigner function, which is a matrix in adjoint color space. Therefore 
it contains many non-Abclian features which are absent from the well known classical 
equation. A notable feature is that, as in the classical case, it contains a term that 
corresponds to the color precession, the non-Abelian analogue to the Larmor precession 
for particles with magnetic moments in a magnetic held. This term was not explicitly 
shown before in the formulation of a Boltzmann equation for the QCD plasma. We find 
that this term is necessary to the gauge covariance of the kinetic equation. 

Acknowledgments 

One of us Q.W. acknowledges financial support of the Alexander von Humboldt- 
Foundation (AvH) and appreciated help from D. Rischke. This work is partially funded 
by DFG, BMBF and GSI. K.R. acknowledges a partial support of the Polish Committee 


Kinetic equation for gluons at the early stage 


6 


for Scientific Research (KBN-2P03B 03018). Stimulating comments and discussions with 
R. Baier, J.-P. Blaizot, E. Iancu and S. Leupold are acknowledged. Our special thanks 
go to X.-N. Wang for his help and interest in this work. 


References 


[1] For a recent review, see eg. S.A. Bass, M. Gyulassy, H. Stocker and W. Greiner; J. Phys. G25, R1 

(1999); H. Satz, Rep. Prog. Phys. 63, 1511 (2000). 

[2] R. Baier, A.H. Mueller, D. Schiff and D.T. Son, Phys. Lett. B502, 51 (2 001). 

[3] For a recent review, see eg. J.-P. Blaizot and E. Iancu, hcp-ph/0101103| , Phys. Rept., to appear. 

[4] L. McLerran and R. Venugopalan, Phys. Rev. D49, 2233 (1994); D49, 3352 (1994); D50, 2225 

(1994); K.J. Eskola, K. Kajantie, P.V. Ruuskanen and K. Tuominen, Nucl. Phys. B570, 379 

( 2000 ). 

[5] S. K. Wong, Nuovo Cimento A65, 689 (1970). 

[6] H.-Th. Elze and U. Heinz, Phys. Rep. 183, 81 (1989); U. Heinz, Ann. Phys. 161, 48 (1985); 168, 

148 (1986). 

[7] H.-Th. Elze, M. Gyulassy and D. Vasak, Nucl. Phys. B276, 706 (1986). 

[8] H.-Th. Elze, Z. Phys. C38, 211 (1988). 

[9] S. Mrowczynski, Phys. Rev. D39, 1940 (1989). 

[10] K. Geiger, Phys. Rev. D56, 2665 (1997). 

[11] M. Gyulassy and X.-N. Wang, Nucl. Phys. B420, 583 (1994); M. Gyulassy, X.-N. Wang and M. 

Pliimer, Phys. Rev. D51, 3436 (1995); X.-F. Guo, and X.-N. Wang, Phys. Rev. Lett. 85, 3591 

( 2000 ). 

[12] R. Baier et al., Nucl. Phys. B483, 291 (1997); B484, 265 (1997); R. Baier et al., Nucl. Phys. 

B531, 403 (1997). 

[13] B.S. DeWitt, Phys. Rev. D162, 1195 and 1239 (1967). 

[14] G. ’t Hooft, Nucl. Phys. B62, 444 (1973). 

[15] L.F. Abbott, Nucl. Phys. B185, 189 (1981); R.B. Sohn, Nucl. Phys. B273, 468 (1986); H. Kluberg- 

Stern and J.B. Zuber, Phys. Rev. D12, 482 (1975). 

[16] H.-Th. Elze, Z. Phys. C47, 647 (1990). 

[17] J.-P. Blaizot, and E. Iancu, Nucl. Phys. B557, 183 (1999). 

[18] J.-P. Blaizot and E. Iancu, Nucl. Phys. B417, 608 (1994); B570, 326 (2000); Phys. Rev. Lett. 70, 

3376 (1993). 

[19] J. Schwinger, J. Math. Phys. 2, 407 (1961); L.V. Keldysh, Sov. Phys. JETP 20, 1018 (1965); G.Z. 

Zhou, Z.B. Su, B.L. Hao and L. Yu, Phys. Rep. 118, 1 (1985); E. Calzetta and B.L. Hu, Phys. 
Rev. D37, 2878 (1988). 

[20] Q. Wang, K. Redlich, H. Stocker and W. Greiner, in preparation. 

[21] T. Altherr and D. Seibert, Phys. Lett. B333, 149 (1994); R. Baier, M. Dirks, K. Redlich and D. 

Schiff, Phys. Rev. D56, 2548 (1977); C. Greiner and S. Leupold, Eur. Phys. J. C8, 517 (1999); 
I. Dadic, Phys. Rev. D59, 125012 (1999); D63, 025011 (2001). 





